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NOTATION 


Coefficient of infinite series, Rauation [41] 
Outer radius of coating 

Effective overheating ratio 

Overheating ratio 

Coefficient of infinite series, Equation [91] 
Wire radius and inner radius of coating 
Cylinder function 

Specific heat at constant pressure + 
Coefficient of infinite series, Equation [68] 
Modified cylinder function 

Group of functions defined in Equation [64] 
Functions defined in Equations [24] and [32] 
Functions defined in Equations [75] and [84] 
Electric current flowing in wire 

Time constant of wire 


Time constant of a constant-current coated wire for a change 
in convective cooling 


Time constant of a constant-current coated wire for a change 
in current input 


Time constant of aconstant-temperature coated wire for a change 
in convective cooling 


Linear combination of cylinder functions, Equations [55] and [56] 
Linear combination of cylinder functions, Equations [33] and [34] 
Function of U from King’s Equation [6] 

Parameter defined in Kquation [21] 

Function defined in Kquation [10] 

Parameter defined in Equation [21] 

Resistivity of wire at temperature T, 

Resistivity of wire at temperature 7, 

Resistivity of wire at temperature 7, 


A radius 


vi 


Function of AP and AT, defined by Equation [18] 
Parameter of the Laplace transform 

Temperature of unheated wire in stream 

Average temperature of the wire 

Reference temperature 

Temperature at radius 7 and time ¢ 

Steady-state temperature at r 

Time 

Fixed parameter 

Flow velocity 


A velocity, characteristic of the flow medium and wire size, 
Equation [6] 


Linear combinations of cylinder functions, Equations [90] and [91] 


Thermal coefficient of resistance 
Parameter defined in Equation [31] 
Parameter defined in Equation [31] 

Ky/ Po Cp, 

Determinant defined in Equations [25] and [76] 
Phase angle 

Group of functions defined in Equation [92] 
Parameter defined in Equation [82] 

3d On cy i 

Phase angle 

Thermal conductivity t+ 


Complex number which replaces s in the inversion integral of 
the Laplace transform 


Strength of heat sources in the wire, Equation [4] 
Parameter defined in Equation [82] 

Density 

Group of functions defined in Equation [77] 


Group of functions defined in Equation [83] 


3 Group of functions defined in Equation [78] 


Go 

db (r,t) Temperature change f 

d (fr, 8) Laplace transform of ¢(r,¢) ft 

u(t) Function defined in Equations [48] and [71] 
@ Circular frequency 

* 


Starred quantities indicate a steady-state value 


} When these quantities have subscripts, 1 refers to wire, 2 to the coating, and 3 to the flow medium. 
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ABSTRACT 


An analysis is made of the unsteady heat flow through a solid cylindri- 
cal wire containing heat sources and a coating in the form of a cylindrical shell. 
The differential equations are linearized and the solutions are found by the 
method of Laplace transforms. The response of a constant-current coated wire 
is obtained for two step-like initiating conditions: (1) a change in convective 
cooling and (2) a change in power or current input. The response of a constant- 
temperature coated wire is obtained for a step-like change in convective cooling. 
An equivalent time constant for each of the three cases is computed. The fre- 
quency responses of the wire to sinusoidal initiating disturbances are derived 
in each case by applying Duhamel’s theorem. Bare wire responses are obtained 


by letting the coating thickness shrink to zero. 


INTRODUCTION 


Hot-wire sensing elements have been used for many years in wind tunnels to measure 
fluctuations in velocity and temperature of the air stream. The wire is usually mounted in 
either aconstant-current or a constant-resistance electronic circuit. In the former case the 
current in the wire is kept constant and the wire temperature and resistance change as the 
rate of convective cooling changes. In the latter case the response consists of a change in 
current which keeps the temperature and resistance of the wire constant. Although there is 
no time lag or distortion in the response of a constant temperature wire it requires a finite 
time to set up temperature changes in a constant-current wire. The distortion in the response 
of a constant-current wire may be corrected by a simple electronic differentiating circuit. 
The elements of this circuit may be adjusted by obtaining a satisfactory frequency response 
to a change in current input. 

Many difficulties arise when attempts are made to use a hot-wire turbulence-sensing 
element in water. One of the most persistent is the formation of a film on the wire, even in 
relatively clean water.! For use in sea water or in other corrosive liquids it would probably 
be necessary to coat the wire with an insulating material. The presence of a film or coating 
not only reduces the wire sensitivity but also increases the time lag in the response, as it 
requires a finite time for temperature changes to penetrate the coating. Therefore, the re- 
sponse of a constant-current coated wire will show more attenuation and distortion than that 
of a bare wire and even the response of a constant-temperature coated wire will show some 
distortion. 


lReferences are listed on page 48. 


The present study was undertaken to determine how the presence of a coating affects 
the wire response and to investigate methods for correcting the response by suitable electronic 
circuits. The responses of a constant-current coated wire are obtained for achange in convec- 
tive cooling and for a change in current input and the results are compared. The response of 
a constant-temperature coated wire is obtained for a change in convective cooling. The corre- 
sponding responses of bare wires are found by letting the coating thickness shrink to zero. 

In all cases the initiating disturbance is a step function and time constants for the step re- 
sponses are obtained. Finally the frequency responses for the three cases are derived from 


the step responses by applying Duhamel’s theorem. 


GENERAL CONSIDERATIONS OF UNSTEADY HEAT FLOW 
IN A COATED HOT WIRE 


In making this investigation, it is necessary to analyze the unsteady heat flow through 
a cylindrical coating and a solid cylinder. The distribution of temperature in the wire and 
coating will depend upon the relative heat capacities of the two media and upon the rate at 
which heat is generated in the wire and dissipated into the flow medium. If the heat conductiv- 
ity of the coating is much smaller than that of the wire there will be an appreciable temperature 
gradient in the coating but only a small temperature gradient in the wire. Furthermore, since 
the rate of cooling depends upon the flow velocity in the immediate vicinity of its surface, an 
angular variation in temperature may also be expected. In the present analysis, however, this 
angular temperature variation has been neglected as well as temperature gradients along the 
wire length. Thus the wire and its coating are assumed to possess cylindrical symmetry and 
the temperature at any point is a function of the radius only. 

Figure 1 shows a sketch of the wire 
and its coating, where a and 6 are, respectively, 
the outer and inner radii of the coating and 6 
is also the wire radius. Let the wire have an 
average resistivity &,, an average temperature 
T,,, and a heating current /. If the wire is 
placed in a uniform flow of velocity U and con- 
stant ambient temperature 7,, a certain steady- 
state temperature distribution will be set up 
within the wire and its coating which will be 
designated as T*(r). It will be assumed that 


the average resistivity in the wire is a linear 


Flow Medium 


function of the wire temperature 


Figure 1 - Sketch of the Wire and Its Coating RE, =R l+e(T, -T,)] J 


where A) is the resistivity of the wire at a reference temperature 7) and a is the thermal co- 


efficient of resistance. The average wire temperature is given by 


9 (b 
T = ai T*(r) rdr [2] 


0 


In the steady state, the temperature in the wire and in the coating are solutions of the 


following differential equations. 
rsb 
Cite 

MBP LED alae |p e 

r dr dr a [3] 
OeSiee Sia 

ve (eae 
r df dr 


where starred quantities refer to steady state values and the subscripts 1 and 2 refer, respec- 


tively, to the wire and coating. The source strength uw? is 


2 
9 ye 
tae 7 b4 K 


[4] 


where x, is the thermal conductivity of the wire. 

At the outer surface of the coating, the wire is cooled by forced convection. The heat 
flux in such cases is proportional to the temperature difference between the coating surface 
and the flow medium. At the center of the wire the temperature must be finite. The two 


boundary conditions are 


dale 
kya ( 2) a {PD [T,*(a) - el =0 
: [5] 
T*() A es 


From King’s study of the heat convection from small cylinders in a moving stream, the propor- 


tionality factor P is the following function of the velocity and properties of the flow medium.? 


p= 2(Y— +1) - $3/ ODS Saale a2 ea) [6] 
27 UY 27 kK, 


where P3rCy > and kK, are, respectively, the density, specific heat, and thermal conductivity 
of the flow medium. In addition to the boundary conditions at r = a and r = O, there are two 
compatibility conditions which apply at r= 6. These conditions require that the temperature 
and heat flux be continuous across the interface. 


[7] 
T,*(b) = T,*(b) = T*(b) 


The solutions to the differential equations which satisfy the boundary and compatibility 
equations are 


* es ilies * 1 Jy (ur) 
T*(r)- 1, +4 = E (8) - Ty 3 en 
[8] 
peers) eels (0) Bye bye Rh ee hoe <= 
2 Ky 0 a B 
where 
Plp*(a)—T1=@ [rec =e 1] [9] 


and J, (u7) is the zeroth order Bessel function of the first kind. The constant Q is defined 
by the relation 


k, #5 J, (ub) = Q Jy (ud) [10] 


If at some time ¢ = 0, an arbitrary time-varying disturbance upsets one or more of the 
parameters U, T., or /, the steady-state conditions will be upset and the temperature will be- 
come a function of the time. The differential equations, boundary and compatibility conditions 


for the unsteady heat flow become 


OT OT 
Le peal: +p? T, (r,t) - 7, + a a eet 
ror or @ n Ot 


Ky @ I + P(t) IT, (a,¢) - T, (¢)] = 0 
i f11] 


OT, 0 T. 
kK, 6 (52) =K,6 ( 1) 9 Wao) 10 (6,2) 
0 ) 
b b 


, 


T,(0,t) 4 = 


T (r,0) = T*(r) 


In these equations 
K K 


a La y =—2 {12] 
PA. 9) iy 


In the foregoing set of equations it will be convenient to express each of the time- 


dependent quantities as the sum of its steady-state and time-dependent increment. Thus 


U(t)=U+AU(t) 1(¢)=1+ AT(¢) 
[13] 


T(¢)=T, +AT_@) EA(r,¢) = T*(r) + h(n, t) 


where U, T., I, etc., are to be considered time-independent parameters unless they are written 
explicitly as functions of time. If the time-dependent increments are assumed to be small com- 


pared with their steady-state values, products and squares of these increments may be neg- 


lected and 
AUC 
2O)= PAN Oe DAL’ maWUAG) 
4a Uy U 
{14] 
2p? 


py? (¢) = wp? + Ap? (t) = pw? + ar AT(t) 


If substitutions are made into Equation [11] and if second-order small quantities are neglected, 
the new set of equations will contain steady-state terms and first-order small quantities. If 
the steady-state equations from Equations [3], [5], and [7] are subtracted, the set of equations 


for the temperature increments ¢, (r,t) and ¢, (7,2) becomes 


it @ [, Cen). 2 2p? ea 1]_1 9%? 
1 Bl Shep fp, (r,t) + So Be) T*(r) 1 Maar 
il at 5) ull db, 
r or OR Wy t 
“20( sel +P, (a,t)=PAT, (t) -(T*(a) - T,) AP(t) [15] 
r 


co) 
n?( $22) = a> ( “1) ; f, (b,t) = , (6,2) 


or b 
f , (0,2) # 


¢ (7,0) = 0 


If the wire is used in a constant-temperature circuit which keeps the average temperature 


constant, there is the additional condition 
b 
i p, (r,t) rdr=0 {16] 
0 


This set of equations is very general and contains all the incremental changes which 
may occur in the operation of a hot wire. The response of a constant-current wire is given by 


@,, (¢), the average temperature change in the wire 
2 b 
be al Br, (Gee) neh 17] 
0 


When the wire responds to a change in convective cooling, A/(¢)=0 and AP and AT. are 
the initiating disturbances. When the wire responds to a change in current input, AP = AT, =0 
and A/ becomes the initiating disturbance. The response of aconstant-temperature hot wire 
to a change in convective cooling is given by the current change A/(¢) and the initiating dis- 
turbances are AP and AT.. For this case Equation [16] applies, as the average temperature 
in the wire does not change. 

In the following sections of this report the problem set up in Equation [15] will be 
solved for three cases in which the initiating disturbances are step functions: (1) a constant- 
current coated wire in which the rate of convective cooling changes, (2) a constant-current 
coated wire in which the current input changes, and (3) a constant-temperature wire in which 
the rate of convective cooling changes. An equivalent time constant and the frequency re- 


sponse will be determined for each case. 


RESPONSE OF A CONSTANT-CURRENT COATED HOT WIRE TO A 
STEP-LIKE CHANGE IN CONVECTIVE COOLING 


When a constant-current coated hot wire is subjected to a step-like change in the rate 
of convective cooling, the wire response is given by the average temperature change in the 
wire. The problem to be solved to obtain the temperature increments ¢, (7,¢) and @, (7,¢) in 


the wire and coating is given by Equation [15] by setting A/ = 0. 


[18] 


0 0 
a,b (2) = «15 ( a fb (b,¢) = b, (0,¢) [18] 


or or 
f,(0,t) 4 & 
(7,0) = 0 


The solutions for the temperature increments as functions of the radius and time are given 
later in Equations [38] and [39]. The wire response to a step-like change in convective cool- 


ing is given in &quation [42]. 
Laplace transforms?'* have been used in obtaining the solutions to the above set of 


equations. The Laplace transform of ¢(r,¢) is defined by the integral 


¢ (7, S) -{ (r,t) eS? dt [19] 
0 


The transformed set of equations is obtained by multiplying each term in Equation [18] by 


eS‘ and integrating over ¢. Then the problem to be solved becomes 


Od pase 
acl (: cy - p? $, (r,8) = 0 
r 


lo a Dama aa Te 
rae a q py (7,8) = 0 


kya (; 2) + Pg, (a,8) -< [20] 


where 


Pp =5 - ,?, q = [21] 
uf) Y 


The solutions to the differential equations are the modified Bessel functions J, and 


Ko: Then 
b, (rs s)=C I, (pr) 


(22] 
bo (7r,s)=A ly (qr) + BK, (q7) 


The constants A, B, and C may be determined by use of the boundary and compatibility con- 


ditions. Finally 


I, (pr) 
Py (r, s) et Hone 
Kk, SA(s) 1, (pd) [23] 
(7, 8) : bD,, (qb,9r) a D,, (qb 
F938) aoe 9 = ? 
Po x, 8A(s) q 10 629997) + Ky oo ‘2 a”)| 
where G is written for 
_ KP 4 | (pd) [24] 
I, (p 5) 
and the determinant A(s) is 
G 
A(s) = qa [a Dy, (92,94) + — Do (9,90) 
297) [25] 
P G 
+ i 76D, (96,94) + sar Doo (9,00) 
2 2 
The functions Diy, (x, y) are the modified cylinder functions, defined by Jaeger* as 
git 
Diy, (es y) = | Up (@) Kg (y) - Kg (2) Ig (y)] [26] 


oz dy 


Some of the properties of these functions and the modified Bessel functions are listed in the 


Appendix. 
The inverse of the Laplace transform is 


Figure 2 - Contour of Integration in the ) Plane 


T +i 0 
AGis—— 
Qt iPofcs 


b(r,r) eA? dd [27] 


The integrand has simple poles at A = 0 and: at 
an infinite number of points A, on the negative 
real axis where A (r,,) = 0. The integration is 
taken along a line, parallel to the imaginary 
axis, which lies to the right of all the singu- 
larities, see Figure 2. The integral from A to 
B plus the integral over a portion C of a circle 
which does not pass through any of the poles 
is equal to 272 times the sum of the residues 
within the contour. As the radius tends to in- 
finity the integral over C tends to zero. Thus 
¢ (r,t) becomes the sum of the residues at the 


poles of the integrand. 


u S Pas (pr) = 1 I, (pr) Ne 
Oe ee (oo can) F ’ Oh aa I, (p 6) i 
my dx 


5 | 78D. (7,97) + £ Dy, (gb, a7) 
Po (r,t) = — ?} 
Do) A (A) se 


co 9 (78, a1) + Dy (6,97) 
is ile pag ech iu Dea EL ING 
n=1 , GA 


FIN A= 


[28] 


At the poleX = 0: p=tp,q=0, and G=Q. The determinant A (\) reduces to 


P= OGe log +) 
(0) aS 3 


Ko 


When g = 0 it may be shown that 
qbDy_ (95,97) + & Do (95,97) =1- @ log = 
Ky Ky b 
At the other poles of the integrand, p=ta,, g=78, and 
A, = nla,” ~ 47) =-yB, 
and G, is defined by 
k,a, 5 J, (a, 5) = G, J, (a, 3) 


It will be convenient to introduce N, (8,7) and N, (8,,r) for the functions 


Ny) (8,7) =7B, 8 DB, 5, ¢B,7) — Dyo %B, 2.tB,7) 
2 


G 
--£ [2,8 Che (6,587) + a Coo (B,, .847)| 


2 


ea | 


G0. Ga WS. Bole. re ¢ Cr (6,2.8,7)| 
LS) 


[29] 


[30] 


[31] 


[32] 


[33] 


[34] 


where (C, 7p (x, y) are the cylinder functions which Jaeger‘ defines as 


gitk 
Ci, (@, y) = [dy (a) ¥ (y) - Yo (w) I, (y) | [35] 
dai ay* 


Some of the properties of these functions are listed in the Appendix. The determinant A (A _) 
may be written in terms of N, (8, a) and N, (8,@) as 


A(,)=-B,2N,(6,2)+— (Ba) =0 [36] 
29) 


Equations [31], [82], and [36] are sufficient to determine the three parameters a,>B,> and 


G,. Finally 
. dA 1 pe BN n2Maea) 
— = = ————_ —> a a 
dr}, MILO) Wee 8 0 Bn 
nh 


G 2 Pan Bs ne AL 
= no Bes PON\\ os Sak Sato a n a,” p2 
eae Ky an Kn 


If substitutions are made in Equation [28] the temperature increments in the two 


J (ur) = J (at) =nla? — pt 
(r,t) = b*(b 0 a OSE 1 
$1 (nt) orn] Bee Dees. Tey [38] 


n=1 


[37] 


regions are 


r Q r oo -yB,t 
by ,t)=¢ mp - oes DAM (6,7) e | [39] 


where the steady-state temperature increment at r = 6, }*(6) is written for 
S PIN = Us @)) Sa PAV 


HO) = —— = [40] 
KA) p_@g ie: Po (1+ 2 toe) 


The coefficients A, are given by 


P 
2 (8, |B - = 1 +E toe $)| 
I SS SS SSS SSS SSS SSS SSS [4 


n 2 PS, 2 
1 OPN MOG. pao oie a pO alt EE (Se, 252 
2 n 0 n n n 
Ky KS Ko a2 Kee 


The wire response is obtained from the average temperature change in the wire. Using 


Equation [17] 
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DN igs 2th a 
rey 22 sro) [1- 34, 224 aie a at 


Ky p? b? n=1 Qa? 
It is reasonable to assume the validity of the following identities: 
by Joonr) _ Jour) 
" Jo (a, 5) Jy (ue) 
n=1 
[43] 


Q 


, ANNG (Bath as 
n=1 
Saari [44] 


If the coating thickness shrinks to zero, the foregoing solution becomes the solution 


for the bare wire. For this case P replaces G,, a = 6, anda, is defined by 


k,a, 6 J, (a, b)=P Jo (a, 5) [45] 


In the coefficient A,, N)(8,,@) = 1 and 


gee 2a,” (P - Q) a 
no a aa ley 20 Lk ee 
ky (a,? ~ p?) G + a, e) 
i 


RESPONSE OF A CONSTANT-CURRENT COATED HOT WIRE 
TO A STEP-LIKE CHANGE IN CURRENT INPUT 


When a constant-current coated hot wire is subjected to a step-like change in current 
input, the wire response is again given by the average temperature change in the wire. The 
problem to be solved to obtain the temperature increments d, (r,t) and ¢, (r,¢) in the wire 
and coating is contained in Equation [15] by setting AP = 0 and AT, = 0. 


als ch oy 2 2 Jo (pr) 1 99, 
r 2( st) + PyOoG) ep a Gd) H Oe 


c . [47] 
od 
x22 ( zt), + Pd (a,t) =0 


[47] 


Ge 


fp, (0,t) # © 


ae 


:) a f, (6,t) = gd, (3,7) 
b 


i) (7,0) =0 


In this set of equations u%*is the following function of the current change A/* which initiates 


the unsteady heat flow. 


2A 1* 


p*= [48] 


* 1 
[? (6) - Al 3] 


The solutions for the temperature increments as functions of the radius and time are given 
later by Equations [59] and [60]. The wire response is given by Equation [62]. 


As in the former problem, the solutions to the above set of equations were obtained by 


the use of Laplace transforms. 


The transformed problem becomes 


@ 
1a (, 9F1 w* J 0 (eT) ¥ 
r ar (- ts). PERO) 7,(ub) 
Od, eared 
ky @ 2) + P ¢,(a,8) = 0 
or 
. [49] 
C) 
“ob fa) Ge fy, $, (0,8) = 6, 0,8) 
f, (0,8) # 00 
where gee =, Pek 
7 y 
The solutions to this set of equations are 
——_ y* yr -Q P 1, (pr) 
ates = b — b 2 
Bs Oo) es eam ae co ao 20) 1, (pb) 
F Jo(ur) — Ig (Pr) [50] 
Jy (ub) 1, (p>) 
( es a Ge aD,, (qr yp (q7,9@) 
$, r,S) = 3 pe + A(s) q 01 qT, ko 00 q759 
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where G is defined in Equation [24] and A(s) is the same determinant which was obtained 
in ®quation [25]. Therefore the poles of this problem are exactly the same as the poles of 
the former problem and the inverse of the Laplace transform is equal to the sum of the 
residues. 

For the pole at the origin p = zp, g= 0, and@=Q. A(0) is defined in Equation [29] 


and the other terms have the values Q? 
2 p2 
2 Kip + Ky 
p2 + pw oe 2 [51] 


y2 Tour) 1 @r) Loa @ Iotur) J (ur) _ ard, (ur) [52] 
p77) Jae), | 1@E)| ¥ 2 UG 2 Jy (ud) 


}2 P a 
qa D,, (gt, 9a)+ — Dy, Gr, ¢a)=— (1+ — log— (53] 
Ko Ko tr 
At the other poles of the integrand where A (A_) = 0; p =ta@, and g=7f8, which were 
determined in the former problem. The determinant A (\,) is defined in Equation [36] as a 


function of N, (8,2), N, (8, a), and P/x, but it may also be defined in terms of new functions 


M, ‘B,, 5), M, (8, 5), and G/x, as 


G 
A(,)=—B, 5 M(B, 5) +—* Mg (B, 6) = 0 [54] 
2 


where 
My (B,7) = ia, 4 Dy (ia, 7, te, 2) + = Doo (itn? ia, 2) 
2 [55] 


=F [250m T Ct, nad + = Coo (a 59| 


iM. (ples = =H 22 eae (a, 7 0, 2) fale Coa, 78, 0)| (56] 
K 
2 


If Equations [54] and [36] are combined, it may be shown that 


My (8B, 5) No (8, 2) = -1 [57] 
If #*(b) is written for 
* (« Bb 2. 2 no lng 7 
p*(6) = - a [58] 
P-@Q se 2 has = 
Ky b 
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the final solutions are 


b, (r,t) = *(d) coke -F [228s anes | 


TG) PIA RED 2 Eo 
[59] 
Se) way ee crwed pads 00 
Za ie Jo (a, TG. O) 
a He log — oo M. ( 2 
[S0) —s7/3 0 
GO ac@) | i. OSE n 60 
#2 Y Law? fon E mame ee 
K b Rite 
2 
The coefficient D, is 
u? 2(G, - Q) A, 


Pao ea? G (A+ ZS) 2 ed) 2) NelPa8 ace D [61] 
Ky tanga te og 


where A, is the coefficient of the former problem, Equation [41]. The coefficients D, may be 
obtained as functions of M, (8,, 5) instead of No (6, 2) by applying Equation [57]. 
The wire response, or the average temperature increment in the wire ¢, (¢), is obtained 


by integrating ¢, (7, ¢) over the cross-sectional area of the wire 


SO aey i 
2Q(1+F) ies 
Ay) ee b*(b le peel (Ae ee 
bu = age #0) y in, = ab [62] 


n=1 


where F is written for 


1 + Eos 2) 


2 
(2 ates [63] 
Q(pu aa log +) “2 


If » 6 is small enough so that the series expansion can be used for Q/k, , F may be written as 


P=9 Late Noes oa _ 5 uf o4 
Fs G-4 tO 225) [64] 


384 
4k, Q+Z2 + — log 2) 


When the heat conductivity of the wire is large compared with that of the flow medium, F is 


very small compared with unity. 
The following identities are assumed to hold 
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1? a 
Sr Jon?) _ Ig (ur) ame 1+ £ log $) Q Jour) pr ae 
Ee es = 
Da Jy (a,b) Jy (ud) i. 


pe IAD) Dog) 
n= 1 kK, (u? b? + — 1h» 2 fey S\N oe 
Ky eee [65] 
es Tee) 14+ flog & 
yi pi eons CEA ates [66] 
a My (B, 5) 1+ log & 


Oy D,=1 [67] 


n= 


For the bare wire G, = P and a, is defined by Equation [45]. Then De becomes 


[68] 
4a? Ti (P - Q)? 


i 2 2 
2 P 

Kip (a, - p?) (« De + =) (or De + = 
Bs] ESS) 


RESPONSE GF A CONSTANT- TEMPERATURE COATED HOT WIRE 
TO A STEP-LIKE CHANGE IN CONVECTIVE COOLING 


If a constant-temperature hot wire is subjected to a step-like change in the rate of con- 
vective cooling, the wire response is given by a current change in the circuit. The problem 


from which 6, (7, fg) po (tr, t), and A/(¢) may be evaluated is obtained from Equations [15] 
and [16]. 


Jy (ur) 
Jo (u 5) 


Ob, 2 2 1 99, 
So Oar) asd sto 7 8e 


1 2 ( So2) 2 fs 
Pr OF or / y 


[69] 


K,@ (2) +P, (a,t) =PAT#* -[T(a)-T,]AP*=S 
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Iba) _ (9% 
iDEN Papa), bh Nota g) Pate) mabe!) 


b,(0,t) # © 


(7,0) = 0 


6 
{ fp, 7,t) rdr = 0 
0 


In this set of equations A/(¢) is contained in y (¢) which is defined as 


2A1(¢) 7 
POs = (710) - ip =z) 


[69] 


[70] 


{71] 


The final solutions for the temperature increments as functions of r and ¢ and for y(¢) are 


given in Equations [88], [89], and [90]. 


Again the solutions of the above set of equations have been obtained by the use of 


Laplace transforms. The transformed problem becomes 
LEG 
r or 


fo) 
281) 92 ,(, hb, (r,s) + p u(s) = Ty Spe 


Jo (ud) 
1a¢ ef) bolts 8) = 0 
Oe A ESE ES 
m0 (2) +P b,(a,8) = 


ag, ad eae hereto 
m0 (St) = nr(S), by (0,8) = 6, (3,8) 


$6, (0,8) # = 


b 
| Gar, s)inar = 0 


where D 


woe Ty 


The solutions to this set of equations for the two regions are 
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[72] 


——_ Ss me ody ie) 2? Ij @e) 
PaO) ars F Jy (ud) H 15 (p8) 


[73] 
S Ko 
aioe 4 aeING) Poy Yo) = 1 Dae enlogere) 
and 
S 
= 74 
oP) s A (s) [74] 
In these expressions 
bT, (pb) 
Fy eset a a Nae (75] 
I, (p 6) 
A(s)=qa [Po (qb, qa) mi a qoD,, (98,90) ] 
[76] 
P 
+ — [Poo a> 9a) ae oD, (gb,90)| 
Ko oO 
and u2 p? 
GO ff 
=o = a [77] 
po +p 


The temperature change as a function of the time is given by the inverse of the Laplace 
transform which is defined in Equation [27]. As in the former problem the inverse is equal to 
the sum of the residues at the poles of the integrand. Again the integrand has simple poles 
at the origin and at an infinite number of points on the negative real axis. 


When A} =0; p=ip,q=0,and H=Q. Then it may be shown that 


ee al (Se es [78] 
% Q 2Q 2k, 
P 
A) =~ (1+ — tog 2+ = [79] 
Ko % 
1 w? Io(ur) — p? Ip (pr) 
p? + 2 E ROD) a a al 
B o oP [80] 


Clee (ee : Tour)! eee ier 
S95 S| = Il) a Ns 
Q 2Q Jo (ub) 25, (ud) 


K 
Doo (qb, qr) -—* 75 D,, (95,97) 
0 


[81] 
aE Slonea eee 
Oh 
At the other poles where A(A,)=0; p= and q=i¢, 
: 
A, =—n(€,2 - w)=-y 62 [82] 
ae 
H 
A eas [83] 
Cn, (ae —p? 
and H, is defined by 
in Go O dh WSO) 2 Tel dh US) [84] 
It will be convenient to introduce the functions Z, (¢,r) and Z, (¢,r) defined as 
K ° . ° 
Gel O) = Don C6,066.0) > = 86 0 Day CEO 8G") 
n 
[85] 


= ae [0 (a b,¢,7) -—2 ¢, O10 (nbtar)| 


Fan eyo tn Mea) 


ul ea, 
pe © | Cay Cn boon?) : nb Oy Cn Balat)] [36] 


The determinant A (A,) may be written as 


AQ,)=-¢,¢2, (¢,4)+ £ Z, (6,4) = 0 [87] 
K 
2 
This expression, together with Equations [82], [83], and [84], is sufficient to determine ¢,, 
(C9 Gas and H_. 
Finally, the temperature increments in the wire and coating are 


Bae) 2) eS er 0 -1) Jour) ur J (ur) 
1+? log 24 P 2102 Q Jy (ud) 2Q Jo (ud) 
b 


OS) % [88] 


+> ats ch Jo (Gr) we Jo (ur) sr Ge HE 
EM Jo (€, 5) @ Jy (ud) 
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S gly etl log 
G Ky b 


by (r,t) = 
eee log 2 Ue 
Ky b or | [89] 


oo Bed 2 
SD Ey hy GaP yon / 


n=1 


The wire response is obtained from the expression 


oo za 2_,2)¢ 
CO) yc aS RUB y7B, ¢ WS ee) [90] 
They LE heya oy Le n=1 
89) % 
where 
A(O 
oes 
(2) 
adn 
r 
n (91] 
-2Z,(¢,@) Yee log 2+ P 
Ko b Fo 
B 7p RES Ce 
P+ 202) 2200 @) eS 
(s R Oa 0,2 Kk, 2 
and 
252 9 2 
Sinai tua gu IG a pin [92] 
H2 H o,, 


The infinite sum in these expressions is an alternating series. The coefficient of the leading 


term B, is positive since Z,) (¢, @) is negative. It is assumed that the following identities 


are valid. 
STB I Semldg GeON eta lu) 
2 £242 AL pea) | O RR 
2 n [93] 
--1|(Ge* 1) Jo (ur) ur J, (ur) 
Q 2Q Jy (ub) 2 Jy (ud) 
y B,Z,(¢,7)=-+- + log * [94] 
Ei Comer: 
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co 2) 327) 
ee | [95] 
— on or Q 2Q 2K, 


The solution for the constant-temperature bare wire may be found by letting a = 0 in the 


foregoing relations. Then 


Mi )eateotis @ [96] 


on 


The characteristic values of €, 6 are determined from the equation 
= Fe = — - [97] 


K 
In the expression for B,, the denominator reduces fo 5, and Z, (¢, a) becomes 


ASA] 
Zy (C,2)=- 2 = 2 [98] 


n 


Then the response of the bare wire becomes 


(Q + 2 pin EF 
oe 3G, 
Qu (t)=-—_S _— 1+ —c 
Niece ate n=1 iy Ee | 
% H, 2H, 2 ky 


Thus the bare|wire responds with a step-function which overshoots its final value. 


[99] 


TIME CONSTANT AND FREQUENCY RESPONSE OF COATED WIRES 


Although the response of a coated hot wire to a step-like change in one of its parameters 
is given by an infinite sum of transient terms, efforts will be made to find out how well a single 
transient term can approximate the sum. An equivalent time constant for the coated wire will 
be determined and its frequency response will be investigated. In addition, sample calcula- 
tions of the time constants and frequency responses will be presented. 

If the response of a system to a step-like disturbance g is characterized by an exponen- 


tial function such as 


f(t)=g fl -e /™) [100] 


the time constant M is defined as the time in which f(¢) attains 1 — e~! of its final value. 
If, on the hand, the initiating disturbance starts with a step but is also a function of 


the time, the response of the system may be found from the solution for the simple step 
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response by applying Duhamel’s theorem.2’” By this theorem the time in g(t) is replaced 
by a fixed parameter ¢’ and the response becomes a function of ¢and¢*% Then f(¢) is ob- 


tained from the integral 
t 
fo 2 | f(t—t,0%) dt’ [101] 
0 
If the response of the system to a sinusoidal disturbance of circular frequency w is required, 
g(t’) may be written as 


g(t’) =g* cos wt’ [102] 


If the values of g(t’) and f(¢,¢’) from Equations [102] and [100] are used in Equation [101], 


the frequency response of the system becomes 


f(t) = g*A(@) cos (wt - &) + cos? 5 et /M [103] 


where the phase angle 6 is 
tan 6 =Mo [104] 
and the amplitude is 


A(@) = cos 6 {105] 


After an appreciable time has passed the transient term becomes very small and the frequency 
response is given by the first term of Equation [103]. When the circular frequency w = 1/M, 


the amplitude has diminished to 0.707 and the response is distorted by a phase angle 6 = 7/4. 


CONSTANT-CURRENT COATED WIRE 


For the case of acoated hot wire an equivalent time constant may also be defined as 
the time in which the wire response attains 1 - e~! of its final value. Thus the time con- 
stant M, of a constant-current coated wire responding to a step-like change in convective 


cooling is, from Equation [42], 


Dy, A, Gag Bn Me xig~1 [106] 
n=1 Qa, 
yB,2M, =1+log A, Ge 
ea [107] 
vlog [1+ 37 Anfant St oe 
n=2 A,G,a,. 
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In most cases the last term is negligibly small. 
The frequency response of the wire may be obtained by writing the initiating disturb- 
ance of Equation [42] as $* cos @t’and by making use of Duhamel’s integral. If the tran- 


sient term is neglected the frequency response becomes 


by, ¢)= 6,7 Alo) cos (wt - 8) [108] 
where a G2 
2, A, Oe 7 Sin, cos 6, 


any S) Sp ae [109] 


n=1 
[110] 
2 
co G, yu 
+ A, 7 sin 6, cos 6, 
iS Qa, 
tan 6, = A JOE ae {111} 


The time constant and frequency response of a constant-current coated wire for a 
change in current input are given by similar expressions in which 4, Gu? /Qe,2 is replaced 


by DG, u?/Qe,? (1+ F). The time constant M, is obtained from the expression 


DC. fe DG een am me 
yByM,=1+ log ed log Lo. So ttl Se eo [112] 
Qe (lua) = DG ice 

Calculations of the time constants and frequency responses for the two initiating con- 
ditions will show how well the response of a constant-current hot wire to a change in current 
approximates the response to a change in convective cooling. In order to make the calcula- 
tions it is necessary to prepare tables of the cylinder functions Cp, (x,y), Co, (4,4), Cy (@¥); 
and C,,(z,y) where y = az/b. These functions are tabulated in Tables 1, 2, and 3 for three 
coating thicknesses: a/b = 1.5, 2, and 3. The functions N, (8, 2), N, (8,7), 4,,D,,M,, and 
M,as functions of P/x, and B, 6 are presented in Tables 4 through 7. For most of these 
calculations the wire parameters chosen were p, Go, DAO ots = 100 x,, and u2b2 = 0.01 
and 0.002. One set of calculations was made for x, = 10 Ko: 


In order to make the set of calculations complete, Table 8 contains similar calculations 
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for the bare wire as functions of 7 (a, — yp”) and P/x,. If P/«, is small, both P/x, and 
Q/x, may be replaced by their series expansions ina, and p. From Equations [45] and [10] 


P a? b2 (: n a,” b2 a’ p4 ) 
— = —_ ——+ oO 00 
Ky 2 8 48 


[113] 


2 2 2 42 444 
ie Jt Giese ee) 
K 8 48 


Then A, and D, from Equations [46] and [68] become the following functions ofa, and pn 


ss det (a, a u?) 6? . (a? + Qi) (a,? Ee u2)b4 


ae) 114 
1 8 96 me 


2 
(a,? - p?) 54 
a ee ee ee, 115 
1 192 (115) 


Thus A, and D, are very nearly unity. Asa, >>a, forn 2 2, the corresponding coefficients 
A, and D,, in the infinite sums are negligibly small. Similarly in the expressions for the wire 


response, Equations [42] and [62], only the leading terms need be considered. In terms of 


a, and » these become 


2 2 2 2) p74 

A Pu Leg Ge Se eres 116] 
1 Qa.? 192 eda 
1 
2 

P yw? a, (a,? = pu?) pe 

Dy eee Soe alee eRe [117] 
Qa, (1+F) 3072 


To a very good approximation the response of a bare wire to a change in convective cooling or 


to a change in current input is given by 


Ra 2_2yt 
a tahoe [118] 


b, t)= 94," [1 - 
To this approximation 


1 Ky b2 
nla? - yp?) 27n(P-@) 
{119] 
4 
mB" 91% =) -T, 


Roal?  pa(b)- Ty += 


where substitutions for 7, P, Q, and p? 62 were made from Equations [12], [9], [10], and [4]. 
This is exactly the result obtained by Dryden and Kuethe® using the assumption that there 
were no temperature gradients in the wire. 

The dimensionless time constants yM,/b? and y M,/b? are plotted against P/x, in 
Figures 3 and 4 and against P/x, in Figure 5. Figures 3 and 4 show the effects of the coat- 
ing thickness and Figure 5 shows the effect of the heat conductivity of the coating. 

In Figure 6 typical curves are drawn for the step responses of a constant-current coated 
wire to a change in current and to a change in convective cooling as functions of t¢/M,. On 
the same graph the response of a bare wire with time constant M, is represented by the dotted 
line. It is clear that the step response to a change in current input is very close to a simple 
exponential function. As the coating thickness decreases and its heat conductivity increases, 
the step response of the coated wire to a change in convective cooling also approaches the 
exponential curve. 

The frequency responses corresponding to these step responses are shown in Figure 7 
as functions of aM). When w = 1/M, the amplitude of the bare-wire response is reduced to 
0.707. The frequency response of the coated wire to a change in current is very nearly the 
same as the bare-wire response but the frequency response to a change in convective cooling 
is lower for all frequencies. The ordinate of the latter curve is approximately 0.707 when 
o=1/M,. 

It is common practice to design a compensating circuit which will correct the distortion 
in the bare-wire response. The elements of this circuit may be adjusted by assuring a satis- 
factory frequency response to a current input since the two time constants are nearly identical. 
As the response of a constant-current coated wire to a step-like change in convective cooling 
is not a simple exponential function a different type of electronic circuit would be needed to 
fully correct the distortion in the response. If the coating is thin, or if Q/x, is only slightly 
smaller than P/x,, the two time constants are nearly equal and the response of the wire to a 
step-like change in convective cooling is well approximated by a simple exponential function. 
In this case the bare-wire compensating circuit can be used in the conventional manner to 
correct the distortion in the wire response. The effectiveness of the compensation decreases 
as the coating thickness increases and as the difference between P and Q increases. 

In hot-wire anemometry the hot-wire element is usually given a particular temperature 
increase over the ambient temperature of the flow. This temperature increment is measured 


in terms of an overheating ratio a, defined as 


Teale ieee Ae 


Wh ey ee ee [120] 
Tie Ty ate a e 
where 
B28, Dl eaGr, =m) [121] 
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Figure 3 - Time Constants for a Constant-Current Coated Wire as Functions of the 
Coating Thickness and P/x, for yu? 6? = 0.01 
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Figure 4 - Time Constants for a Constant-Current Coated Wire as Functions of the 
Coating Thickness and P/x, for y? 6? = 0.002 


(Kk, = 100 Ky Py Py = Po “Py 


A | 


AGE 


ee eS 


0.6 


0.06 0.08 0.1 


0.03 0.04 


= 0.01) 


= 2, ye 0? 
27 


Effect of the Thermal Conductivity of the Coating Material on the 
(a/b 


Time Constants of Constant-Current Coated Wires 


Figure 5 - 
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Figure 6 - Responses of a Constant-Current Coated Wire to a Step-Like Change 
in Current and Convective Cooling 
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Figure 7 -- Frequency Responses of a Constant-Current Coated Wire 
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which is the average resistivity of the unheated wire. If use is made of the relations given 
in Equations [8] and [9] and if the approximation is made that the average temperature of the 
wire is the same as its surface temperature and Q = 4% Ky 2 b?, then a,, is the following 


function of P and Q. 


Q E A Le leet) 
kK b 
— [122] 
IPa=@ (is Eh log 2) 
Ky b 


Values of a,, are included in Tables 4 through 8. 

Although the overheating ratio of a coated wire may be reasonably high, the tempera- 
ture at the outer surface of the coating may not be very much above the ambient temperature 
of the flow. In this case the wire will be more sensitive to temperature fluctuations than to 
velocity fluctuations. If an effective overheating ratio a, is defined in terms of the surface 
temperature of the coating rather than the wire, it may be shown that 


Ee EO es ee ea [123] 


Pomp. fe LE a 
1 Uy ee OE 
2 
Although a, cannot be measured directly it is a more significant parameter than a, for a 
coated wire. 


CONSTANT-TEMPERATURE COATED WIRE 


In a similar manner the step and frequency responses of ‘a constant-temperature coated 
wire may be found. From Equation [90] the wire response to a step-like change in convective 


cooling may be written as 


oo we 2 
BOD TS he yBn "| [124] 


n=1 


where yw* contains the initiating function and other constants. The equivalent time constant 


M, is defined by the equation 


= B, -y(82-B,2)M 
FORE SR osetia BS ni a eM [125] 


The corresponding frequency response is 


u(t) = p* A(@) cos iwt — $) [126] 
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where Pe 
> B, sin@, cos 6, 
enn a ae ee [127] 
1-5 B, sin20, 
1 


n= 


z 0° 2 co 2 
Ato) = (: - > B, sin? ) a S B, sin@, cos *) {128] 


n=1 n=1 


tan @, = —® [129] 


In preparing to make the calculations for the time constants it was observed that 
'$*(b) was very small if k, 2100 «, and if u2b2<0.01. Also k,/o,, is small for the first 
few values of n. Therefore, most of the calculations were made for Ky = Kp/ On =0% 
(x,/x,) 5, = 0 and $*(5) = 0. For this approximation the problem is independent of yu? b? 
and 


Mau AA LIC, (CB Gin) 
2 [130] 


Z(C,r)= 2 Cy, Cy 8.057) 
The functions C'),(z,y) and Cy, (z, y) are tabulated in Tables 1, 2, and 3 for three coating 
thicknesses: a/b =1.5, 2, and 3. The data for obtaining the time constants are presented 
in Tables 9 through 11. Curves showing how the time constant depends upon the coating 
thickness and P/k, are presented in Figure 8. ' 

If the thermal conductivity of the wire is not large compared with that of the coating, 
k,/o, and (x,/k,) 5, cannot be neglected. In order to find out how the thermal conductivity 
of the coating affects the wire response, calculations of the time constant were made for 
a/b = 2, u? 62 = 0.01, and for kK, =100 ky and 10 Kae The data for these calculations are 
given in Tables 12 and 13. The results for x, = 100 x, were almost the same as the results 
for xk, = ©. The curves of Figure 9 show how the time constant depends upon P/x, and 
K,/ Ko. 

The response of the bare constant-temperature hot wire to a step-like change in con- 
vective cooling is itself a step-function with an initial overshoot and appears to be quite 
unlike the response of the coated wire. Inthis case a time constant has no meaning and the 
frequency response is flat. If Q/x, is replaced by its series expansion in Equation [99], 
the bare-wire response is 
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Figure 8 - Time Constants for Constant-[emperature Coated Wires as Functions of the 
Coating Thickness and P/k, (x, = ~) 
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Figure 9 - Effect of the Thermal Conductivity of the Coating on the Time 
Constants of Constant-Temperature Coated Wires 
(a/b = 2, p2 % = 0,01) 


33 
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Sear Va + i 5 
a3 2, 128 nn x 


w(t)=w* {14+ {131] 


If P/«, is small compared with unity, the values of H,,/x, are only slightly less than P/x, 
and the roots €, 6 of Equation [97] are very nearly the zeroes of J, (€, 5). Therefore 
Ky ge b2/2 H,, even for n = 1, is very large compared with the other terms of the denominator 


and Equation [131] is given approximately as 


OP — ome, -H7)t 
p(t)=w*l1 ne aa Y CL SD ERE [132] 
Ky = Be b2 
Initially 
0) =-S ay fi+ 
WO ex ( + ne) [133] 


Therefore a constant-temperature bare hot wire responds to a step-like change in convective 
cooling with a step-like response which has a small overshoot of magnitude P/4«,, which 


decays rapidly. The time constant for the decay of the overshoot is obtained from the 


equation 
wo —n(€7—"2)M 
S eat Ga (134] 
2 72 
n=1 or b 
and Fi 


1M ~ 0.0834 


The response of a constant-temperature hot wire is proportional to the heat flux at the 
wire surface 7 = 5. If the wire has no coating there is no time lag in the response (except for 
the small overshoot) and the frequency response is flat. If the wire has a coating, a step- 
like change in the heat flux at the outer surface is no longer a step function when the flux 
reaches the inner surface and there is an attenuation and a phase distortion inthe wire re- 
sponse. This becomes apparent if the step and frequency responses of a coated wire are 
examined for particular values of the wire parameters. Such responses are shown in Figures 
10 and 11. As the infinite sums converge rather slowly for small values of ¢ and for large 
values of w, there is an uncertainty of a few percent inthe step response for small values of 
t and in the frequency response for aM, >1. As the coating becomes thinner and P/x, 
becomes larger the step response becomes steeper, less like an exponential function and 


more like a stepfunction. As long as there is a thin coating there will be some delay in the 
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Figure 10 - Response of a Constant-Temperature Coated Wire to a Step-Like 
Change in Convective Cooling 


(a/b= 2, Ky = 0% P/ky = 9.62) 
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Figure 11 - Frequency Response of a Constant-Temperature Coated Wire 


(a/b=2, Ky = Pky = 9,62) 
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response and a meaningful time constant. When a = 5 there is no longer a time delay and the 
response becomes a step. 

This interpretation is borne out by the frequency response which is more characteristic 
of a step response than of an exponential response to a step-like initiating function. Unlike 
the frequency response of a constant-current coated wire which approaches the dotted line 
on Figure 11 from below, the frequency response of a constant-temperature coated wire is al- 
ways above this line and approaches the line A (w) = 1 as the coating shrinks to zero. 

These results show that the type of electronic circuit used to correct the response of 
a constant-current hot wire would not be effective in the case of a constant temperature coated 
wire. Although it might be possible to design a circuit for particular values of the wire param- 


eters, it would be difficult to adjust the components of the circuit for different flow conditions. 


SUMMARY AND CONCLUDING REMARKS 


In setting up the problem for the unsteady heat flow in coated hot wires the assumption 
was made that the initiating disturbance as well as changes in the wire parameters were small 
compared with their steady-state values. Second order small quantities were neglected and 
the differential equations were then linear inthe temperature change and other incremental 
quantities. Solutions were obtained for step-like changes in the initiating conditions and 
equivalent time constants were determined for the wire response. From the step responses 
of the wire the frequency responses were found by applying Duhamel’s theorem. 

Although the time constants for the step response of a constant-current bare wire for 
a change in convective cooling and for a change in current input are nearly identical, they 
become different when the wire has a coating. Whereas a single exponential function is a 
good approximation for the step response of a coated wire to a change in current input, it 
becomes a poor approximation for the step response to a change in convective cooling, par- 
ticularly for a large coating thickness. This is seen in Figure 6 where the two responses 
are compared. The time constant for the change in convective cooling M, is larger than that 
for a change in current M;. The wire also has a better frequency response for changes in 
current input than for changes in convective cooling, as shown in Figure 7. 

When the coating is thin and the two time constants are nearly equal the wire response 
may be corrected by a simple differentiating electronic circuit. The compensation becomes 
less effective as the coating thickness increases, particularly if the elements are adjusted 
by correcting the response to a current input. Therefore, another method of setting the ele- 
ments in the compensation circuit should be found. For perfect compensation a more compli- 
cated electronic circuit must be used. 

The response of a constant-temperature hot wire is quite different from that of a 
constant-current wire. The constant-temperature bare wire responds with a step to a step-like 
change in convective cooling and it has a flat frequency response. When the wire has a coat- 


ing there is a delay in the response and the step function is replaced by a more gradual rise, 
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see Figure 10. The freauency response of a coated wire falls off with increasing frequency 
but the curve is always higher than the corresponding frequency response curve for the 
constant-current wire, see Figure 11. 

As the step response of a constant-temperature coated wire is quite different from a 
simple exponential function the simple differentiating circuit used to compensate a constant- 
current wire would be wholly inadequate. From the data given here it would be possible to 
design a suitable circuit but it might be difficult to determine a method of setting the circuit 


elements for different operating conditions. 
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TABLE 1 


Values of the Functions Cp 9(#,y), Coy (#,Y), ®Cy o(# ¥), and £C,,(z,y) 
for a/b = 1.5 


[> [es [ean rn a 


0.254262 0.671056 | —0.611646 | 0.156798 
0.247477 0.364740 | —0.568144 | 0.254431 
0.237502 0.221317 | —0.505004 0.342003 
0.224575 0.133442 | —0.424674 0.416057 
0.209002 0.072114 | —0.330257 0.473669 
0.191148 0.026366 | —0.225410 | 0.512567 
0.171430 | —0.008802 | —0.114168 | 0.531225 
0.150300 | —0.035969 | -0.000823 | 0.528921 
0.128237 | —0.056625 | +0.110280 | 0.506769 
0.105732 | -0.071713 0.240180 | 0.462697 
0.083273 | —0.081907 


0.061337 | —0.087747 
0.040369 | —0.089715 
0.020781 | -0.088271 


—0.05524 +0.001514 
— 0.05387 0.006835 
—0.05202 0.01188 
0.04974 0.01660 
—0.04706 0.02098 
— 0.04403 0.02499 
—0.04067 0.02861 
— 0.03705 0.03182 
— 0.03319 0.03460 
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TABLE 2 


Values of the Functions C), (2, y); Co (2, y), C,9(#,y), and xC,, (z, y) 
for a/b = 2 


Cooley) | Coley) | 2C, (ey) } 2C,(z,y) 


0.438311 | +1.551182 
0.429507 | +0.715881 
0.415065 | +0.412752 
0.395332 | +0.244683 
0.370779 | +0.132801 
0.341984 | +0.051205 
0.309624 | -0.010953 
0.274449 | -0.058947 
0.237262 | -0.095648 
0.198903 | —0.122799 
0.160372 | —0.141577 
0.122034 | —0.152895 
0.085152 | —0.157538 
0.050312 | —0.156252 
0.018175 | —0.149780 
—0.01069 —0.13888 
—0.03581 —0.12433 
—0.05685 —0.10693 
—0.07357 —0.08747 
—0.08586 —0.06675 
—0.09376 —0.04553 
—0.09741 —0.02453 
—0.09704 —0.00442 
— 0.09301 +0.01419 
— 0.08576 +0.03081 
—0.07579 +0.045016 | -0.17831 | -0.40825 
— 0.06366 +0.056491 | —0.25847 | -0.36012 
— 0.04994 +0.065028 | -0.32788 | -0.297538 
—0.035255 0.070529 | 0.38382 | -0.22313 
~0.020187 0.072996 | -—0.42402 | -0.13986 
—0.005309 0.072540 | -—0.446840 | —0.050995 
+0.008849 0.069355 | —0.451456 | +0.039802 
0.021811 0.06372 | —0.437719 | 0.128918 
0.033174 0.05598 | -0.406225 | 0.212786 
0.04261 0.04655 | -0358274 | 0.288036 
0.04987 0.03585 | —0.295855 | 0.351691 
0.05481 0.02435 | -0.221497 | 0.401184 
0.05737 0.01252 | —0.138229 | 0.434560 


0.05758 | +0.000813 
0.05683 | -0.00486 
0.05554 | —0.01034 
0.05373 | -0.01559 
0.05144 | -0.02056 
0.04870 | -0.02521 
0.04554 | -0.02950 
0.04201 | -0.03340 
0.03815 | -0.03689 
0.03399 | -0,03993 
0.02960 | -0.04252 
0.02501 | -0.04463 
0.02028 | -0.04625 


—0.04090 +0.001341 
—0.04029 0.00536 
— 0.03929 0.00925 
—0.03792 0.01299 
—0.03619 0.01652 
— 0.03415 0.01984 
—0.03179 0.02290 
—0.02917 0.02569 
—0.02630 0.02817 
—0.02321 0.03033 
—0.01995 0.03217 
— 0.01655 0.03365 


—0.626379 | 0.094872 
—0.596045 | 0.186047 
—0.546769 | 0.269975 
—0.480416 | 0.343391 
—0.399498 | 0.403445 
—0.307070 | 0.447815 
—0.206619 | 0.474795 
—0.101917 | 0.483366 
+0.00310 0.47286 
0.10456 0.44484 
0.198658 | 0.399365 
0.281971 | 0.338614 
0.351479 | 0.265031 
0.404718 | 0.181546 
0.43988 | +0.091479 
0.45585 | -0.001597 
0.45230 | —0.09400 
0.42962 | —0.18208 
0.38899 | —0.26235 
0.33226 =| -0.33172 
0.26192 | -0.38737 
0.18095 | -0.42724 
0.09274 | -0.44974 
+0.00095 | -0.45403 
—0.09060 | -0.44000 
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for a/b = 3 


(u2b? = 0,01 and 0.002, 


0. ae 0.43845 
0.32005 
0.15221 


0.54040 
0.41941 
0.24781 
0.03674 


0.694655 
0.680540 
0.657395 
0.625774 
0.586434 
0.540309 
0.488487 
0.432176 
0.372678 
0.311316 
0.249477 
0.188489 
0.129632 
0.074094 
0.022936 


—0.15756 
—0.15052 
— 0.13833 
~—0.12181 
— 0.10185 
— 0.07941 


0.09380 
0.09026 
0.08340 
0.073 64 
0.06146 


0.61389 
0.62008 
0.67076 


0.34483 
0.47364 
0.58167 
0.64656 


2.059820 
0.937865 
0.525879 
0.294562 
0.138926 
+0.024580 
—0.062813 
-0.130196 
— 0.181308 


—0.627412 
—0.600130 
—0.555772 
—0.495965 
—0.422898 
~0.339237 
—0.248028 
—0.152572 
— 0.056309 


0.218516 | 


— 0.243325 
—0.257009 
— 0.260689 
—0.255458 
—0.242428 


— 0.00126 
+0.02847 
0.05483 
0.07718 
0.09506 
0.10814 


—0.000855 
— 0.01889 
—0.03533 
~ — 0.04965 
—0.06139 


K, = 100 k, 


Pe Oe 
2’ "1 py 


TABLE 3 — Values of the Functions C,)(z,y), Co,(@,¥), @C,9(@,y), and «C,, (x,y) 


eC, ( (x,y) 


0.084347 
0.165507 
0.240416 
0.306252 
0.360550 
0.401328 
0.426976 
0.437121 
0.430183 


TABLE 4 — Time Constants M, and M, for a Constant-Current Coated Wire a/b = 1.5 


my) 
PoP. 


TABLE 5 


Time Constants M , and M, for a Constant-Current Coated Wire a/b = 2 


Da a a 
(uw? b* = 0.01 and 0.002, kK, = 100 kp, Pip, = Pa ep, 
0.39819 0.90242 0.90651 25.247 


0.34659 
0.26403 


0.60183 
0.55114 
0.15557 0.54445 | 5.59943 1.7688 
0.09363 0.54204 | 10.4202 1.4225 
0.02779 0.53658 | 38.6138 1.1639 


lanl fees mle PP LL 


0.01 -—0.69791 | —0.18251 | 1.25526 | —0.0387 
: : —0.68502 | —0.35204 | 2.67234 : 
: : —0.62961 | —0.51693 | 4.59778 b 
. : -0.53190 | —0.66686 | 7.52235 ; 
: ( 0.39491 | —0.79327 | 12.8559 
—0.22382 | —0.88728 | 26.957 


25.2805 | 24.6910 
6.5238 | 5.9861 
3.0387 | 2.5766 
1.4269 
1.1288 
0.9207 
0.7693 


1.38916 
2.50494 


6.5083 
3.0174 


0.57376 
0.51870 
0.43047 
0.31420 
0.24763 
0.17670 
0.10247 


0.28108 
0.31496 
0.38563 
0.44627 
0.46823 
0.48325 
0.49063 


0.19595 
0.48577 
1.0750 
2.2725 
3.4035 
5.4698 
10.5339 


—0.12041 
—0.28801 
—0.45299 
—0.60545 
—0.73572 
— 0.83600 


— 0.64834 
— 0.65041 
—0.60904 
—0.52445 
—0.39934 
—0.23859 


0.89144 | —0.0941 
2.3026) | -0.1937 
4.16521 | ~0.2558 


6.92667 | - 0.2941 


} 
11.7906 | ~6.3161 
23.827 | -oans 
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TABLE 6 


Time Constants M, and M, for a Constant-Current Coated Wire a/b = 2 


Dh 
(u°b* = 0.01, kK, = 10 Ky, p1¢ 2 °p,) 


Ea oe EN NC EE 


0.07012} 0.59564 0.20365 0.13676 | 1.0112 24.655 | 0.669 
0.13042} 0.54003 0.27942 0.41393 | 1.0442 | 0. : 9.9678] 0.184 
0.23133 | 0.45075 0.36546 0.97293 | 1.0976 | 0. : 2.9691} 0.094 


0.37347) 0.33277 0.43477 2.0904 | 1. : : 1.4220} 0.062 
0.46025] 0.26504 0.45988 SoS} |) le : : 1.1273} 0.053 
0.55770} 0.19271 0.47751 4.9556 | 1. : c 0.9224) 0.047 
0.66596} 0.11685 0.48704 9.1698 | 1. : d 0.7738} 0.042 
0.78520) 0.03854 0.48802 | 30.3881 | 1. ; : 0.6614} 0.038) 


TABLE 7 


Time Constants M, and M, for a Constant-Current Coated Wire a/b =3 


(2 b2 = 0.01 and 0.002, k 1 = 100K 


Q? Pye Py = Py? Pp 


100.840 
25.840 
11.918 
6.9818 | 5.8631 | 1.992 | 0.347 


4.6220 | 3.7365 Ean | 
r 

100.899 | 97.456 2.634 | 2.235 
25.890 | 23.447 0.726 | 0.536 
11.984 
7.0939 
4.7956 
3.5058 
2.6626 


1.2228 
1.5956 
2.4090 
4.3225 
12.0243 


1.12579 
0.65376 
0.52739 
0.47733 
0.44751 


0.27619 
0.24584 
0.19703 
0.13252 
0.05583 


0.50563 
0.52068 
0.54574 
0.58084 
0.62598 


0.10503 
0.12004 
0.14506 
0.18008 
0.22513 
0.28020 
0.34530 


0.30069 
0.27809 
0.31670 
0.35930 
0.39183 
0.40822 
0.40529 


0.16270 
0.32184 
0.61908 
1.1249 
2.0205 
3.9117 
10.7255 


0.55445 
0.51844 
0.46041 
0.38328 
0.29087 
0.18784 
0.07935 
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TABLE 8 


Time Constants M, and M, for a Constant-Current Bare Wire a/b = 1.0 


(2 b? = 0.01 and 0,002) 


Cae NO oa WE 


0.007213 | 0.999448 69.4443 
0.009824 | 0.998796 | 1. 51.0203 
0.012841 | 0.998043 | 0.999999 | 39.0624 
0.016266 | 0.997298 | 0.999997 | 30.8640 


0.020101 | 0.996231} 0.999995} 24.9998 | 


0.0018008 | 0.999800 | 1. 277.178 | 277.778 
0.0032026 | 0.999575 | 1. 156.250 | 156.250 
0.0050063 | 0.998999 | 1. 100. 100. 
0.0072130 | 0.998448 | 0.999999 ) 69.4444} 69.4444 
0.0098240 | 0.997796 | 0.999998) 51.0203} 51.0204 


TABLE 9 
Time Constant M, for a Constant-Temperature Coated Wire a/b = 1.5 


(Ky =o, Py “py =Po “Dp? 


yM,/0? 


0.23100} 3.0 | 9.39 | 1. 383 | 0.13566 
1.2205 | 3.4 | 9.52 | 1. 505 | 0.11083 


2.9169 : ; : 605 | 0.09326 
4.2730 : ; : : 0.08025 
6.7868 c : ¢ : 0.07022 
10.7293 : ; : ¢ 0.06214 
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TABLE 10 


Time Constant M, for a Constant-Temperature Coated Wire a/b = 2 


(Ky =") 91°, =poe, ) 
1 1 Py 2 Po 


0.09905 
0.6873 
1.4513 
2.4695 
3.8843 
6.0139 
9.6204 
17.3917 
49.445 


TABLE 11 
Time Constant ,, for a Constant-Temperature Coated Wire a/b = 3 


(Ky =o, PA, =P Day 


0.27003 
0.72302 


1.3136 
2.1057 
3.2186 
4.9087 
7.8429 
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TABLE 12 
Time Constant M, for a Constant-Temperature Coated Wire a/b = 2 


2,2 
(K, =100k,, p,¢, =p,¢, , pb“ =0.01) 
1 OPEN a Ba? 


0.0025041 | —0.31014 | —0.012142 | 0.10962 
0.0025048 | -0.27470 | —0.060158 | 0.70078 


0.0025055 | —0.23725 | —0.096833 | 1.46931 
0.0025063 | -0.19864 | —0.12391 | 2.4952 
0.0025072 | -0.15987 | —0.14258 | 3.9240 
0.0025081 | -0.12133 | -0.15374 | 6.0825 
<6 | 0.0025091 | —0.084270 | —0.15820 | 9.7621 


TABLE 13 


Time Constant a for a Constant-Temperature Coated Wire a/b =2 


kK, = 10K), p,¢ Py = Pye Po rp 22 = 0.01) 


0.025228 | -0.31484 | —0.022931 | 0.20394 
0.025292 | —0.27703 | —0.71172 | 0.82213 
0.025366| —0.23718 | -0.10764 | 1.6338 
0.025449 | -0.19624 | -—0.13412 | 2.7338 
0.025541} -0.15530 | -0.15178 | 4.3002 
0.025644 | —0.11480 | -0.16158 | 6.7557 


0.025756) -—0.07610 | —0.16436 | 11.2313 


0.09960 | 0.007970 
0.09303 | 0.02683 
0.08321 | 0.04319 
0.07072 | 0.05662 
0.05623 0.06684 
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APPENDIX 
PROPERTIES OF BESSEL FUNCTIONS AND CYLINDER FUNCTIONS 


The Bessel differential equation may be written as 


=e) n? 

1d + 

es sel Nt pire ws -0 
Lee ? ae 


If the sign of ¢ is positive the solutions are the nth order Bessel functions of the first and 
second kinds J, (x) and Y (x). If the sign of ¢ is negative the solutions are the modified 
Bessel functions /, (x) and K, (x). If mis an integer the two sets of functions are related by 


the equations 


1, (iz) = Fi)” J, (@) 
K, (fiz) = Gar** ZU, (@) F 2 ¥, (@)] 


Some of the important properties used in this report are listed below.2"® As J, (a) and Y, (x) 
have the same properties, Z, (x) will be used here to represent either function or any linear 


combination of these functions. 


Z, (@) = Z,_,(@) hs Z, (t)=-Z,,,(@) + = Z,, (a) 


1, (a) = 1,_, (@) re I, (e) = 1,,,(@) + = I, (x) 
Ky (@)=-K,_,(@)+4 K, @)=-K,,,(2) + 2 K, () 
Zy (#) =-Z,(@) = Z_,(@) 


Iy(z)= 1,(@)= I, (2) 


Kf (@) = - K, (@) = - K_, () 
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The following integral relations are needed 


[2 Zo(erae=2 2,00) 


| Z, (x) de =- Z, (z) 


ar Z, (az) Zo (82) - Bx Z, (Br) Z (az) 


|. Z, (ax) Z,(B x) dx = ra 
a safe 


[- be! (az) dz = _ [Z” (az) + Z}(az)] 


When the boundary conditions are applied the solutions of some of the differential equa- 
tions of this report are the cylinder functions Cik (a, y) and the modified cylinder functions 


Dix (a, y). Jaeger* defines these functions as 
Coo (ts ¥) = Iq (a) Yq (y) - Yq (z) Jy (y) 


j+k 
C., (2, y) = : C, (a, 
pW alae? OO y) 


Doo (ty) = Ig (t) Kg (y) - Kg (2) I (y) 


j+k 


Dix (t, y) = el ay Doo (ty) 


These two sets of functions are connected by the relation 
gitk D jy, ia, ty) = - 5 Ci, (ty) 


The modified cylinder functions have the following properties 
1 


Dj krg\tsy) = Dix (x,y) - ; Di ry 2) 


iL 


Diy (x, z) = =D (x, 2) = : 
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Similar properties for the cylinder functions are 


Cita, bry) = — Cip(ery) a Ci41,46%Y) 


Ci pyg(tsy) = - Cy (ay) - - Op nace) 


2 

C 9 =-C 9 =— — 

10 (% 2) o1 (% 2) or 
C.. (x,y) C,, (2, y) - C,, (a, y) C (ay) ee 
00 MY) Oy MP 01 S79 Y) O49 UY mle) 


In the limit as 2 and y approach zero 


Dy (2a, ty) = oo Coo (ty) = - log # 
ix D, 9 («, ty) = ar tC (2,y) =1 
y Do, (tz, vy) = ins y Coy (t,y)=-1 


vy Diy (tx, cy) = a vy Cy (z, y) =0 
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